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It is argued that small values of the neutrino masses may be due to anthropic selection eflFects. 
If this is the case, then the combined mass of the three neutrino species is expected to be ~ 1 eV, 
neutrinos causing a non-negligible suppression of galaxy formation. 



I. INTRODUCTION 

The major ingredients of the Universe are dark energy, 
f2A ~ 0.7, and non-relativistic matter, ~ 0.3. The 
latter consists of non-baryonic dark matter, Qo ^ 0.25, 
baryons, ~ 0.05, and massive neutrinos, il^ > 0.001. 
The fact that is comparable to Qm is deeply puz- 
zling; this is the notorious coincidence problem that has 
been much discussed in the recent literature. The only 
plausible explanation that has so far been suggested is 
that r^A is a stochastic variable and that the coincidence 
is due to anthropic selection effects. Anthropic bounds 
on the cosmological constant derived in [1-4] were fol- 
lowed by anthropic predictions [5-8] suggesting values 
not far from the presently observed dark energy density. 
Although controversial, such anthropic arguments have 
been bolstered by the discovery of mechanisms that may 
be capable of creating ensembles with different parame- 
ter values in the context of both cosmic inflation [9-11] 
and string theory [12-15], and have been applied to other 
physical parameters as well [16-33]. 

Perhaps equally puzzling are the "coincidences" Qo ^ 
fls and Qb ^ ^i^- These three matter components are 
relics of apparently unrelated processes in the early Uni- 
verse, and it is very surprising that their mass densities 
are comparable to one another. The mass density of neu- 
trinos is fli, = (mi//94eV)/i~^, where is the combined 
mass of all three neutrino flavors. In this paper, we will 
investigate the possibility that rrii, is a stochastic variable 
taking different values in different parts of the Universe 
and that the observed value is anthropically selected. 

Before delving into details, let us briefly outline the ar- 
gument. A small increase of can have a large effect on 
galaxy formation. Neutrinos stream out of the potential 
wells created by cold dark matter and baryons, slowing 
the growth of density fluctuations. As a result, there 
will be fewer galaxies (and therefore fewer observers) in 
regions with larger values of m^. If the suppression of 
galaxy formation becomes important for nii, ^ rhi,, say, 
then values ^ will be rarely observed because the 
density of galaxies in the corresponding regions is very 
low. Moreover, unless the underlying particle-physics 
model strongly skews the neutrino mass distribution to- 
wards values near zero, values rrii, <^ to^ are also unlikely 
to be observed, simply because the corresponding range 
of mi/-values is very small. A typical observer thus ex- 



pects to find ^ m^, i.e., a mild but non- negligible 
suppression of galaxy formation by neutrinos. 

II. PROBABILITY DISTRIBUTION FOR 

To make the analysis quantitative, we define the prob- 
ability distribution V(m^)dm,y as being proportional to 
the number of observers in the Universe who will mea- 
sure rui, in the interval dnii,. This distribution can be 
represented as a product [5] 

V{my) = V*{my)nobs{rni,). (1) 

Here, V*{mi,)dm,y is the prior distribution, which is pro- 
portional to the comoving volume of those parts of the 
Universe where rrii, takes values in the interval drriu, and 
nobs(jni,) is the number of observers that evolve per unit 
comoving volume with a given value of rriu- The distri- 
bution (1) gives the probability that a randomly selected 
observer is located in a region where the sum of the three 
neutrino masses is in the interval dm^. 

Of course we have no idea how to estimate riobs, but 
what comes to the rescue is the fact that the value of 
does not directly affect the physics and chemistry 
of life. As a rough approximation, we therefore assume 
that nobs{mu) is simply proportional to the fraction of all 
baryons that form stars, which we approximate by the 
fraction -Fm(wj^) of all matter that collapses into galaxy- 
scale haloes (with mass greater than M = 10^"^ AIq), 

nobsirrii,) cx FM{rn^). (2) 

The idea is that there is some average number of stars 
per unit mass in a galaxy and some average number of 
observers per star. The choice of the halo mass scale is 
based on the empirical fact that most stars are observed 
to be in large halos. 

The prior distribution Vt:{m^) depends on the exten- 
sion of the particle physics model which allows neu- 
trino masses to vary and perhaps on stochastic processes 
during inflation which randomize these variable masses. 
Some candidate prior distributions will be discussed in 
Section III. 

The fraction of collapsed matter FM{'my ) can be ap- 
proximated using the standard Press-Schechter formal- 
ism [34]. We assume a Gaussian density fluctuation 
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field (5(x, t) with a variance cr(M, t) on the galactic scale 

(M = IQI^Mq), 



P((5, cx exp 



(3) 



A collapsed halo is assumed to form when the linearized 
density contrast 6 exceeds a critical value 6c determined 
by the spherical collapse model. As detailed in Appendix 
A, this corresponds to dc ~ 1-69 around the present epoch 
and Sc « 1-63 in the infinite future [4]. Using the Press- 
Schechter approximation, we obtain 



/>oo 

FM{t) oc P[6 >5c]= / P{6, t)d6 = erfc 



V2aM{t) 



(4) 

The collapsed fraction Fm thus grows over time as the 
rms density fluctuations a increase. 

Let us now quantify the effect of neutrino masses on 
this process. For the small scale M that we are consider- 
ing, assuming a flat Universe, this fluctuation growth is 
well approximated by 



aM{x) 



1 + 2^a(/.)Ga(x) 



ctm(O), (5) 



as shown in Appendix A. The functions ^a, Ga and p 
are defined below. Here we have replaced t by a new time 
variable 



Pa 



(1 + z)^n„ 
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i.e., the dark-energy-to matter density ratio — we will 
consider several values of x below, corresponding to the 
infinite future {x = oo), the present epoch {x = 7/3, 
our default value) and redshift unity (a; = 7/24). The 
function 



Ga{x) 
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where a = 0.795 and 
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1.43728, 
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describes how in the absence of massive neutrinos, fluc- 
tuations grow as the cosmic scale factor a as long as dark 
energy is negligible (Ga (.'/;) « .x^/^ oc a oc (1 -I- z)~^ for 
a; -C 1) and then asymptote to a constant value as t ^ oo 
and dark energy dominates {G\{x) G^o as x ^ oo). 

We are (;onsidering the case where r/v,^ varies from place 
to place whereas the physics that determined the amount 
of baryons and cold dark matter per photon is the same 
everywhere, so the neutrino fraction is given by 



f = ^ — 

i im Phc + Pv 



1 + ^ 
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where pbc denotes the non- neutrino density, i.e., that of 
baryons and cold dark matter, and mbc ~ (4.75±0.30)eV 
gives the measured amount of such matter per neutrino. 
In other words, increasing the neutrino mass from zero 
will increase the total matter density per photon by a 
factor Pm/Pbc = (1 - /i/)~^- 



-1/3 



1/3 



-f (1+^eq) (10) 



is the factor by which the Universe has expanded between 
matter-radiation equality at a; = x^q (when fluctuations 
effectively start to grow) and dark cmergy domination at 
a: = 1 (when fluctuations gradually stop growing). Since 
massive neutrinos boost the matter density by a factor 
(1 — /y)~^, they delay vacuum domination until the scale 
factor is larger by a factor (1 + fuY^^ and also, in the 
approximation that neutrinos are nonrelativistic at the 
matter-radiation equality epoch (valid for rriu ^ 1 eV), 
cause matter-radiation equality to occur occurs earlier, 
when the scale factor is smaller by a factor [1 + f^). We 
thus have 



Aa{U)^{1-U)-^'^A^{Q). 



(11) 



Finally, neutrinos with nonzero mass suppress the galaxy 
density through the exponent p{fv) in equation (5), 
which is given by [35] 



V25 - 24/, - 1 



(1-/.) 
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(12) 



and drops from unity for = to smaller values as 
increases. 

In summary, equation (5) shows that the galaxy fluc- 
tuation evolution (Jm{x) depends on the three cosmo- 
logical parameters Aa(0), (7m(0) and To study the 
galaxy density as a function of neutrino fraction fi, us- 
ing equation (4), we thus need to measure ^a(O) and 
o'm(O) from observational data without making any as- 
sumptions about /i/. We cannot do this using the val- 
ues of and Zeq reported by, say, the WMAP team 
[36], since these assume that /j, = in our part of the 
Universe; if > here, then matter-radiation equal- 
ity occurred earlier. We therefore repeat the Monte 
Carlo Markov Chain analysis reported in column 5 of 
Table 3 in [37], measuring Aa(0) and aM(0) from the 
WMAP microwave background power spectrum [38] com- 
bined with the Sloan Digital Sky Survey (SDSS) galaxy 
power spectrum [39]. These measurements arc indepen- 
dent of fu since this is a free parameter in the analysis 
and therefore effectively marginalized over. This gives 
Aa(0) = 3057 ± 502, aM(0) « 0.000579 ± 0.000064^. 



^ As our galactic scale, wc take M = 10^^ Mq, specifically a top- 
hat smoothing scale of _R = l.Sfe^^Mpc. This corresponds to 
length scales about 100 times smaller than the matter-radiation 
equality scale where the matter power spectrum turns over. 
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The above-mentioned rribc-value was measured using this 
same MCMC analysis. We will use the central values 
for our main analysis and quantify the effect of the un- 
certainties in the discussion section. Equation (5) thus 
shows that for f,y = 0, fluctuations grow by a factor 
1 4- 1.5Aa{0)Ga{x) w 4700 by the present epoch, which 
we take to be a; = OA/ilm « 0.7/0.3 ~ 2.3, giving 
ctm ~ 2.7. In the infinite future x — > oo, fluctuations 
will have grown by a factor 6600, giving (Jm ~ 3.8. The 
basic reason that neutrinos have such a dramatic effect 
is that these growth factors are so large, implying that 
even a modest change in the exponent p{fi,) makes a large 
difference. Taylor expanding equation (5) in f,^ gives 

fTM(a;,/.)«aM(x,0)e-^(-)-^'' (13) 

for U < 1, where j{x) = 0.61n[l -I- 1.5^a(0)Ga(x)] - 
4/3 « 3.7 for the present epoch and 7 ~ 3.9 for the infi- 
nite future. Although equation (13) is quite a crude ap- 
proximation, underestimating the suppression, it shows 
that small changes in A\ or x are unimportant since they 
affect this exponential fluctuation suppression only loga- 
rithmically. 

The effect of neutrino free streaming on the galactic 
density is illustrated in Fig. 1 (top), which shows that 
the suppression is non-negligible already for rrii, ~ 1 eV. 
We use equation (5) in our calculations for the plots — 
the approximation equation (13) was merely to provide 
qualitative intuition for the effect. 

The probability distribution V{mi,) is shown in Fig. 1 
(bottom) for power law priors 

■P4m^)(xm[!, (14) 

with n ranging from O to 4. For n > 1, these distributions 
are peaked at ^ 2 eV, while in the case of a flat prior, 
7^*(toi/) = const, the expected values are ~ 1 eV. 
This is also seen in Fig. 2, where the distribution for a 
flat prior is shown using a logarithmic scale for m^. 

In this discussion, we have assumed that /,y ^ 1, that 
is, nil, ^ 10 eV. Very heavy neutrinos (with 3> 1 
MeV) would annihilate well before nucleosynthesis and 
cause no problems with structure formation. If all neu- 
trinos were heavy, neutrons would be stable, leading to an 
equal number of protons and neutrons. As a result, most 
of the matter would end up in helium instead of hydrogen. 
This lack of hydrogen would clearly suppress riobsi'm-u) for 
observers like us who rely on long-lived (hydrogen burn- 
ing) stars and water-based chemistry. Moreover, heavy 
neutrinos would not be able to blow off the envelope in 
supernova explosions. This means that heavy elements 
formed in stellar interiors would not be dispersed to form 
planets and observers. The possibility of the electron 
neutrino being light and one or two others very heavy 
is allowed anthropically, but it is already ruled out by 
the neutrino oscillation experiments, which constrain the 
mass differences to be within 0.05 eV. 

For m„ > 100 eV (but < 1 MeV), keeping aU 
other physical parameters fixed, neutrinos would have 
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FIG. 1: The upper panel shows the factor by which the neutrino 
fraction (dashed curve) suppresses the current fluctuation ampli- 
tude cruj (upper solid curve) and consequently the galaxy number 
density na (lower solid curve). The lower panel shows the resulting 
probability distribution for the neutrino mass sum for priors 
with n = 0, 1, 2, 3 and 4, peaking from left to right, respectively. 



sufficiently low thermal velocities to act approximately 
as cold dark matter, thereby allowing galaxy-size halos 
to form. However, the baryon fraction in these halos 
would be strongly diluted, and it is therefore far from 
clear whether they would be able to cool and form self- 
gravitating baryonic disks, let alone stars or observers, 
with an efficiency comparable to that in our observable 
universe. In other words, the calculation of anthropic 
constraints on very large neutrino masses becomes essen- 
tially equivalent to the calculation of an anthropic upper 
bound on the dark matter abundance. We will not at- 
tempt to address this issue here, but simply assume that 
the number of observers noi,s{my) is strongly suppressed 
for my > 10 eV. 

We have also assumed that there are N^j — 3 stable 
neutrinos. Generalizing our result to iV^ > 3 is straight- 
forward: as long as the masses are low enough for neu- 
trino infall to be negligible, the galaxy number density 
depends only on the total neutrino mass density, which 
for standard neutrino freezeout is proportional to the 
sum of the neutrino masses. If the neutrinos are un- 
stable on cosmological timescales, they suppress fluctua- 
tion growth only before decaying, with their decay radi- 
ation redshifting away to gravitationally negligible levels 
within a few expansion times. 
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FIG. 2: The probability distribution for the neutrino mass sum for 
flat prior (n = 0). The dark/red tails contain 5% probability each. 
The dotted line shows the lower limit 0.05 eV from atmospheric 
neutrino oscillations [41-44]. 



III. PRIOR DISTRIBUTION 

Following [45] , wc shall now discuss possible modifica- 
tions of the standard particle physics model that could 
make neutrino masses variable. For early work on how 
masses of elementary particles can vary randomly in the 
context of stochastic gauge theories, see [46-48]. 

Dirac-type neutrino masses can be generated if the 
Standard Model neutrinos mix through the Higgs dou- 
blet VEV $ to some gauge-singlet fermions v^, 

go.^^v'-yt (15) 

The couplings gafi will generally be variable in string the- 
ory inspired models involving antisymmetric form fields 
Fa interacting with branes. (Here, the index a labels dif- 
ferent form fields.) changes its value by AF^ — 
across a brane, where qa is the brane charge. In the low- 
energy effective theory, the Yukawa couplings gaj3 become 
functions of the form fields, 

ffo/3 = 5S+E5S,^ + -, (16) 

where the summation is over all form fields, the coeffi- 
cients g^^^ , g'^p are assumed to be numbers ~ 1 , and Mp 
is the effective cutoff scale, which we assume to be the 
Planck mass. 

In such models, closed brane bubbles nucleate and ex- 
pand during inflation [49], creating exponentially large 
regions with different values of the neutrino masses. 
When Fa changes in increments of (/q, changes in 
increments of 

Am, ~ ^qalMl. (17) 



To be able to account for neutrino masses ^ 1 eV, we 
have to require that Am, < 1 eV, that is, 

qa < 10-"M2, (18) 

for at least some of the brane charges. Such small val- 
ues of the charges can be achieved using the mechanisms 
discussed in [12, 13, 50, 51]. 

It should be noted that the Higgs potential and the 
Higgs expectation value $ will also generally depend on 
Fa- Moreover, each field Fa contributes a term F^/2 to 
the vacuum energy density pA , and regions with different 
values of Fa will generally have different values of pA- 
However, in the presence of several form fields with suf- 
ficiently small charges, variations of all these parameters 
are not necessarily correlated, and here we shall assume 
that there is enough form fields to allow independent vari- 
ation of the relevant parameters. We can then consider 
a sub-ensemble of regions with m, variable and all other 
parameters fixed. The probability distribution for to, 
that we calculated in Section II corresponds to such a 
sub-ensemble. 

Let us now turn to the prior distribution V*{m^). The 
natural range of variation of Fa in Eq. (16) is the Planck 
scale, and the corresponding range of the neutrino masses 
is < rriv^ ^ (Here, the index i labels the three 
neutrino mass matrix eigenvalues.) Only a small fraction 
of this range corresponds to values of anthropic interest, 
'T^i/ ^ 10 eV. In this narrow anthropic range, we expect 
that the probability distribution for Fa after inflation is 
nearly flat [52], ^ 

dV* = const • dFidF2..., (19) 

and that the the functions gapiFa) are well approximated 
by linear functions. If all three neutrino masses vary 
independently, this implies that 

dV^ = const • dm';^Umf> dm'i;'' . (20) 

The probability for the combined mass m, = ^ toI'-* to 
be between to, and m, -I- dm, is then proportional to the 
volume of the triangular slab of thickness ~ dm, in the 
3-dimensional mass space, 

dV^ cx ni^drav (21) 



^ A very different model for the prior distribution was consid- 
ered by Rubakov and Shaposhnikov [40]. They assumed that 
T'priori^) is a sharply peaked function with a peak outside the 
anthropic range A and argued that the observed value of X 
should then be very close to the boundary of A. We note that 
this is unlikely to be the case for the neutrino mass, since it is ob- 
served to be comfortably inside the anthropically allowed range. 
If the model of [40] applied, the peak of the full distribution 
would most likely be in a life-hostile environment, where both 
Vpriori^) and riabsi^) a^re very small. In the case of the neu- 
trino mass, this would mean that the number density of galaxies 
is very low. This is not the case in our observable region, indi- 
cating that the model of [40] does not apply. 
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Alternatively, the neutrino masses can be related to one 
another, for example, by a spontaneously broken family 
symmetry. If all three masses are proportional to a single 
variable mass parameter, then we expect 

dVt. oc dm^. (22) 

Let us now assess how well the predictions derived from 
the prior distributions (21) and (22) agree with observa- 
tions. We first consider the distribution (21), correspond- 
ing to independently varying neutrino masses. The most 
probable value of m,y for n — 2 is ^ 3 eV, and we 
expect both the neutrino masses and mass differences 
to be '--^ 1 eV. This expectation, however, is in conflict 
both with neutrino oscillation experiments suggesting 
Am^ < 0.05 eV [41-44] and with astrophysical bounds 
[36, 37] which indicate a combined mass rrii, ^ 1 eV. 

For a flat prior distribution (22), the most probable 
value is nii, ^ 1 cV. If nii, is close to this value, then the 
three neutrino masses must be nearly degenerate, with 
Am^ <C m^. This could be interpreted as a sign of a 
family symmetry. A 90% confidence level prediction for 
nil, based on this distribution can be obtained as outlined 
in Section II. This gives 

0.1 eV < < 5 eV. (23) 

The lower bound in (23) is somewhat stronger than the 
bound from the neutrino oscillation data [41-44] (mi, > 
0.05 eV), while the upper bound is somewhat weaker 
than the current astrophysical bounds {e.g., [36, 37, 53- 
56]). Note that the strength of current astrophysical 
bounds is limited not by statistical errors but by sys- 
tematic uncertainties in non-CMB data. For instance, 
the recently claimed evidence for m,^ > [57] may result 
from underestimated galaxy cluster modeling uncertain- 
ties. 

We finally mention the possibility that the right- 
handed neutrinos have a large Majorana mass ^ 
$. In this case, small neutrino masses can be generated 
through the see-saw mechanism, 

~ g^<i>yMR. (24) 

If A/fl is variable, say, within a range Mfj ^ Mp, then 
its most probable values are likely to be ~ Mp, and the 
prior distribution will be peaked at very small values of 
m„ ~ 10-6 eV. 

This discussion suggests that the most promising sce- 
nario with variable neutrino masses is the one with Dirac- 
type masses determined by a single variable mass param- 
eter. It yields a flat prior distribution for m,^, Eq. (22), 
and the prediction (23) at 90% confidence level. 

After we submitted the original version of this paper, 
Jaume Garriga pointed out to us that see-saw-type mod- 
els can yield cosmologically interesting prior distributions 
for nil, if the Majorana mass is restricted to the range 
Mr < Mjj'""^-' < Mp. Assuming first that Mr is a fixed 
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FIG. 3: Same as Figure 2, but showing the robustness of the results 
to changing various assumptions. We have changed the baseline cal- 
culation from Figure 2 (heavy black curve) by evaluating the galaxy 
density in the infinite future x — > oo (dotted red/grey curve) and at 
redshift unity, x = 7/3(1 -I- z)^ = 7/24 (dot-long-dashed blue/grey 
curve), decreasing the density threshold to 5c = 1-5 (short-dashed 
magenta/grey curve), lowering the primordial fluctuation ampli- 
tude on the galactic scale by 25% (solid blue/grey curve), including 
the baryon correction as per Eisenstein & Hu (long-dashed green 
curve) and additionally including the neutrino infall correction as 
per Eisenstein & Hu (dot-dashed cyan/grey curve). 

constant, while g is variable, Eq. (24) yields the distribu- 
tion 

(X dg (X m ' dm. (25) 

This would give a somewhat smaller predicted neutrino 
mass than the distribution (22) that we used in most of 
our calculations. 

The distribution (25) applies up to nimax ^ /Mr. 
In order to have nimax ^ 0.1 eV, we need Mr < 10^'^ 
GeV. 

If both Mr and g are variable, then, assuming a flat 
prior for Mr, Eq. (25) still applies, but now nimax ~ 
^2/^j{rnax)^ SO wc need Mjj"""^ < 1013 GeV. An attrac- 
five feature of this scenario is that the increment of nii, in 
Eq.(17) gets suppressed by an additional factor ^/Mr, 
and Eq.(18) gets replaced by a much weaker constraint 
qa/M^ < (Mh/1013 GeV)-i. 



IV. DISCUSSION 

In conclusion, we have found that the small values of 
the neutrino masses may be due to anthropic selection. 
If so, then the most promising model appears to be the 
one with a flat prior distribution, 'P*(to;^) = const. The 
range of m-i, predicted in this model, Eq. (23), has in- 
teresting implications for both particle physics and cos- 
mology. On the particle physics side, neutrino masses in 



6 



this range arc nearly degenerate, suggesting extensions 
of the Standard Model involving a spontaneously broken 
family symmetry. On the cosmological side, a combined 
neutrino mass of > 1 eV has a non-negligible effect on 
galaxy formation. This means that it must be taken into 
account in precision tests of inflation that measure the 
shape of the primordial power spectrum by combining 
microwave background and large-scale structure data. 

Let us close by discussing the importance of the as- 
sumptions we have made and outlining some open prob- 
lems for future work. The purpose of this brief paper and 
the prediction of equation (23) is merely to demonstrate 
that anthropic selection effects may be able to explain 
the neutrino masses, and much work needs to be done to 
place this hypothesis on a flrmer footing. 



A. Robustness to approximations and 
measurement errors 

To quantify the robustness of our results. Figure 3 
shows how the probability distribution for rrii, changes 
when various assumptions are altered. 

First of all, the parameters A\ and (Tm(0) that we used 
have non-negligible measurement uncertainties. We see 
that lowering crM(O) by 25% (by about twice its mea- 
surement uncertainty) lowers the m,^-prediction slightly. 
Changing Aa by within its observational uncertainty has 
an even weaker effect since, as we saw, it enters only 
logarithmically. Altering the galactic scale M affects gm 
and hence the results only weakly, because of the flatness 
of the dimensionless power spectrum k^P{k) on galactic 
scales. 

Second, our calculations involved various approxima- 
tions. We used the Press-Schechter approximation with 
density threshold 6c ~ 1.69 as per Appendix A; lowering 
this to account for post-virialization infall as discussed in 
[7] is seen to raise the mj^-prediction slightly. Our fluctu- 
ation growth treatment of equation (5) is highly accurate 
in the limit of small mass scales M and low baryon frac- 
tion Qb/flm -C 1, agreeing with a CMBfast [58] numer- 
ical calculation to within a few percent. Figure 3 shows 
that for the observed baryon fraction ilb/rim ~ 0.15, 
switching to an exact treatment of baryon effects makes 
virtually no difference. The cosmic expansion eventu- 
ally slows neutrinos enough for them to start clustering 
on galaxy scales, and if this happens before dark energy 
domination, then it reduces 7 equation (13). Since a 
small fraction of the neutrinos will be in the low tail 
of their velocity distribution, there is a slight infall cor- 
rection even for the low m^-range we have considered, 
and Figure 3 shows that this increases our m,^-prediction 
slightly. Finally, we have used the cutoff value of x « 7/3, 
which amounts to using the reference class of observers 
in galaxies that have formed by now. Figure 3 shows 
that if we ask instead what would be observed from a 
random galaxy among all galaxies that ever form (set- 
ting X = 00), then the mj^-prediction increases slightly. 



Conversely, it also shows that considering only observers 
in galaxies that formed by redshift unity decreases the 
TOy-prediction. In conclusion, Figure 3 shows that al- 
though many of these assumptions make marginal dif- 
ferences, none of them affect the qualitative conclusions, 
since they all shift the predicted probability distribution 
by much less than one standard deviation. 



B. Effects of other parameters 

The standard models of cosmology and particle physics 
involve of order 10 and 28 free parameters, respectively. 
In order to apply anthropic constraints to them, it is 
crucial to know both which of them can vary, and what 
the interdependencies or correlations between them are. 
It is likely that at least some of the cosmological pa- 
rameters (the baryon-to-photon ratio, say, via baryogc- 
nesis) are determined by particle physics parameters in 
a way that we have yet to understand, and many par- 
ticle physics parameters may in turn be determined by 
a smaller number of parameters or vacuum expectation 
values of some deeper underlying theory. A proper anal- 
ysis of anthropic predictions should therefore be done in 
the multi-dimensional space spanned by all fundamental 
variable parameters. 

Such correlations between parameters must ultimately 
be taken into account not only for computing the theo- 
retical prior Vt: of equation (1), but also when computing 
the factor Uobs in this equation, which incorporates the 
observational selection effect. The reason is that strong 
degeneracies are present which can in many cases offset a 
detrimental change in one parameter by changes in oth- 
ers. For instance, suppressed galaxy formation caused 
by increased to^ can to some extent be compensated by 
decreasing p\, by increasing the dark- matter-to-photon 
ratio or by increasing the CMB fluctuation amplitude Q 
above the value ~ 10~^ that we observe [24, 32] — if any 
of these three parameters can vary, that is. In the present 
paper, we have merely considered the simple case where 
all relevant parameters except rrii, {i.e., the comoving 
densities of baryons and dark matter, the physical den- 
sity of dark energy, and the fluctuation amplitude Q) arc 
kept flxed at their observed values, with no account for 
scatter due to variation across an ensemble or from mea- 
surement uncertainties. A more detailed study of this 
issue is given in [59] and shows that our present results 
for are rather robust to assumptions about pA- 

This is closely related to the issue of how much infor- 
mation one wishes to include in the factor Uoba in equa- 
tion (1) [32, 60]. One extreme is including only the ex- 
istence of observers, the other extreme is including all 
available knowledge (even, say, experimental constraints 
on m,y). As one includes more such information, the an- 
thropic factor becomes progressively less important, and 
the calculation acquires the flavor of a prediction rather 
than an explanation. In the context of a multiparameter 
analysis, the question is whether to use the measured val- 
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ues of other parameters (in our case non-neutrino param- 
eters) or marginalize over them. Our fixing non-neutrino 
parameters at their observed values is therefore equiva- 
lent to factoring in the information from the measure- 
ments of these parameters. 

Arguably the most interesting outstanding question is 
whether the fundamental equations that govern our Uni- 
verse do or do not allow physical quantities such as the 
neutrino masses to vary from place to place. Calcula- 
tions of anthropic selection effects may prove useful for 
shedding light on this. In any case, for quantities that 
do vary, the inclusion of anthropic selection effects such 
as the one we have evaluated in this paper is clearly not 
optional when calculating what the theory predicts that 
we should observe. 
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APPENDIX A: GROWTH OF LINEAR 
PERTURBATIONS 

In this Appendix, we derive and test the approximation 
given by equation (5) for how small-scale matter fluctua- 
tions grow in the presence of radiation, cold dark matter 
and neutrinos. There are two reasons why this simple 
approximation complements an exact "black box" calcu- 
lation with CMBfast [58] or a nearly exact approximation 
with the Eisenstein & Hu fitting software [61]. First, for a 
qualitative argument like the one we make in this paper, 
it is desirable to have a simple intuitive understanding 
of the underlying physics that includes only those com- 
plications that really matter for the argument. Second, 
neither CMBfast nor the Eisenstein & Hu package were 
designed to be valid for extreme cosmological parame- 
ters such as those corresponding to the infinite future, 
and indeed break down in this limit. 



The ACDM case 
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FIG. 4: The LCDM growth function Ga is shown as a function 
of cosmic time, first growing as the scale factor a oc x^^^, then 
asymptoting to 1.44 as dark energy halts the fluctuation growth at 
a: > 1. The middle panel shows the exact result of equation (Al) 
(solid curve), our approximation given by equation (7) (green thick 
dashed curve), the Carroll, Press & Turner approximation [62] (red 
dotted curve) and the power law approximation Ga = f^m^^ [63] 
(blue long-dashed curve). In the bottom panel, the various ap- 
proximations have been divided by the exact result, showing that 
equation (5) is accurate to better than 1.5% for all x. The top 
panel shows the collapse density threshold (5c (x) dropping from 
1.6865 early on to 1.62978 in the infinite future. 



X oo (when Ga Goo)- Figure 4 shows that this 
approximation greatly improves on the standard Carroll, 
Press & Turner [62] and power law fits for our present 
purposes, since these were designed to be accurate only 
in the past and present and have the wrong limiting be- 
havior in the future as a; ^ oo, and IIa — > 1. 
For flat models, we have fim = l/(l + a;), ilx = x/{l + x) 
and X = (r^m^ ~ so in terms of the standard linear 
growth factor D = Gk/x^^^, the three approximations 
shown in Figure 4 are 



D 



-1/37 



(A2) 



For a flat Universe with only pressureless matter (dark 
and baryonic) and a cosmological constant, the growth of 
density fluctuations is given by (5 oc Gt^{x)^ where [7, 64] 



Ga{x] 



dy 



(Al) 



D « 

2 



350rJ„i 



UOnT + (209 - f^m)^^m + 2 



and 



We find that our Ga{x) fit defined by equation (7) is 
accurate to better than 1.5% for all x and becomes ex- 
act both in the limits a; — > (when Ga x^^^) and respectively. 



-0.21 



1 -1 



(A3) 



(A4) 
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2. Including radiation 

Early on, dark energy was negligible but radiation was 
gravitationally important, causing density fluctuations to 
grow as 5 oc G^(x), where [65] 

o / \ 1/3 
G,{x) = l + -( — ) . (A5) 

Since both Gj{x) and G\{x) accurately describe the 
growth during the matter-dominated epoch x^q <^ x <^ 
1, with Gy cx G\ oc x^^"^ during this period, we can com- 
bine them to obtain the approximation 

G{x) « 1 + ^AaGa{x), (A6) 

which is accurate for all x. Here the constant A\ is de- 
fined by equation (10). In essence, fluctuations grow as 
(5 oc a oc x^/^ between matter domination {x = Xgq) and 
dark energy domination {x = 1), giving a net growth 

1 /3 

of Aj\ = Xcq ■ Equation (A6) shows that they grow 
by an extra factor of 1.5 by starting slightly before mat- 
ter domination and by an extra factor Ga(oo) « 1.44 by 
continuing to grow slightly after dark energy domination. 



4. The collapse density threshold 5c 

In the top panel of Figure 4, we have numerically com- 
puted the collapse density threshold Sc as a function of 
cosmic time x, defined as the linear perturbation the- 
ory overdensity that a top hat fluctuation would have 
had at the time when it collapses. We see that it varies 
only very weakly with time (note the expanded vertical 
scale in the figure), dropping from the familiar cold dark 
matter value 6^0) = (3/20)(127r)2/3 « 1.68647 early on 
to the limit Sc{oo) = (9/5)2-2/3(7^ ~ 1.62978 [4] in the 
infinite future. This calculation neglects the effect of neu- 
trinos. Since their effect is to contribute a gravitationally 
inert component just as dark energy, we will ignore their 
effect on Sc{x). assuming that they merely cause a slight 
horizontal stretching of the curve (which is seen to be 
almost constant anyway). 



3. Including neutrinos 

As shown by [35], the result 5 oc a is generalized to 
^ oc a^' when a fraction of the matter is clustering in- 
ert and remains spatially uniform. The new exponent 
p < 1 is given by equation (12) where, in the case that 
we are focusing on here, the inert fraction is the neu- 
trino fraction f^.^ This motivates our approximation in 
equation (5), which simply generalizes equation (A6) by 
introducing the neutrino-dependent exponent p{fi,). 

We have tested this approximation by comparing equa- 
tion (5) with exact results using the CMBfast software 
[58] and the semianalytic approximation of Eisenstein fc 
Hu [61], finding excellent agreement (to within a few per- 
cent) with both in the small-scale limit for a; ~ 1 and 
negligible baryon fraction. In the distant future limit 
X ^ oo, both CMBfast and the Eisenstein & Hu software 
break down, since they were not designed to be accurate 
for such unusual parameter values (Oa = 1, etc.). For 
the parameter ranges of interest to us, there are small 
corrections for the effects of both baryons and neutrino 
infall, which we quantified in Figure 3 in the discussion 
section. 



^ The result is more general [61, 66], and applies also when the 
inert density components correspond to dark energy or spatial 
curvature. If we let Qm denote the density fraction that is not 

inert (that clusters), then the approximation to equation (12) 

3/5 

given by p Q^i is quite accurate, being exact both for f2m = 
and to first order in (1 — f2m) for all 1 — Qm 1. This is the 
familiar result that dlnS/d\na ~ f^m® ■ 
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